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Emptiness formation probability in the domain-wall 
six-vertex model 

F. Colomo and A.G. Pronko 



Abstract. The emptiness formation probability in the six-vertex model with 
domain wall boundary conditions is considered. This correlation function al- 
lows one to address the problem of limit shapes in the model. We apply 
the quantum inverse scattering method to calculate the emptiness formation 
probability for the inhomogeneous model. For the homogeneous model, the 
result is given both in terms of certain determinant and as a multiple integral 
representation. 



1. Introduction 

A special kind of fixed boundary conditions, the so-called domain wall boundary 
conditions, was first considered for the six-vertex model by Korepin in seminal 
paper [1]. In [2], Izergin showed that the partition function of the model on the 
finite lattice can be found exactly in terms of certain determinant; see also paper [3] 
for details. It was later shown in [4,5], by studying the thermodynamic limit, that 
the free energy per site is different with respect to the case of periodic boundary 
conditions. This fact hints at spatial separation of phases (e.g., ferroelectric order 
and disorder), which is confirmed both numerically [6-8] and analytically [9,10]. 

To get some details about the phase separation phenomena, e.g., to find the 
shape of the spatial curve separating the phases, or limit shape, one has to know 
some appropriate correlation function. The problem of computing the correlation 
functions in the 'domain- wall' six- vertex model has been addressed in papers [11-14] 
where some correlation functions near the boundary were found. 

In this paper, we continue the study of correlation functions of the six-vertex 
model with domain wall boundary conditions. Specifically, we consider here a par- 
ticular non-local correlation function, the emptiness formation probability (EFP). 
This function describes the probability of having a set of consecutive horizontal 
edges along a given column, all in a given state; we consider here the case when the 
set starts from the top boundary and extends inside the lattice. This correlation 
function allows one to address the problem of limit shapes in the model [15]. 

To compute EFP, we follow the lines of papers [12, 14, 16] where the quantum 
inverse scattering method (QISM) [17, 18] and some facts from the theory of or- 
thogonal polynomials were used. Mostly following ideas (as well as notations) of 
[14], we represent EFP in certain determinantal form, which is shown here to be 
also equivalent to some multiple integral. This last representation recalls analogous 
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multiple integral representations for correlation functions of quantum spin chains 
[19-22]. 

The paper is organized as follows. In the next Section we start with giving some 
definitions and fixing some notations. The quantum inverse scattering method in 
application to the model is considered in Section 3. The core calculation of EFP 
for the inhomogeneous model is contained in Section 4. The homogenous limit 
is performed in Section 5, where the main result is given given both in terms of 
certain determinant and as a multiple integral representation. Section 6 is devoted 
to discussion of equivalent multiple integral representations for EFP. 

2. Some definitions and notations 

2.1. The model. The six-vertex model is a statistical mechanics model in 
which the local states are associated with edges of a square lattice, and the Boltz- 
mann weights are assigned to its vertices. The states can take two values, which 
are often denoted by arrows pointing along the edge. Among the sixteen possible 
arrow configurations around a vertex only six are allowed (having nonzero Boltz- 
mann weights), with equal number of incoming and outgoing arrows. In this paper 
we consider the model on a lattice having both N rows and N columns ('the N x N 
lattice') with the boundary states fixed in a special way: all arrows on the left 
and right boundaries are outgoing while on the top and bottom boundaries all ar- 
rows are incoming. Such a model is called the six-vertex model with domain wall 
boundary conditions. 

In the six- vertex model with invariance under reversal of all arrows there are 
three possible values for Boltzmann weights at each vertex, usually denoted as a, b, 
and c. To use the quantum inverse scattering method (QISM) for calculations we 
will consider the inhomogeneous version of the model, in which the weights of the 
vertex being at the intersection of fc-th horizontal line and a-th vertical line are 

Oq/c = a(Ac(, I/fc), bak = b{\a,Vk), Cak = Ci (2.1) 

where 

a(A, i^) = sin(A — + 77), h{\,v) = sn\{\ — v — r}), c = sin(2r/) (2.2) 

and we enumerate vertical lines (labelled by Greek indices) from right to left, and 
horizontal lines (labelled by Latin indices) from top to bottom. The parameters 
Ai, . . . , Aat are assumed to be all different; the same is assumed about i^'s. After 
applying QISM we set these parameters equal within each set, Aq. = A and V]^ = v 
and, without losing generality, we can assume that z/ = 0. In this way we obtain 
the homogenous model; the indicated procedure will be referred to as homogeneous 
limit. 

The partition function of the inhomogeneous model is defined in a standard 

way as the sum over all possible configurations, each configuration being assigned 
its Boltzmann weight, whcih is the product of all vertex weights over the lattice, 

JV JV 

C a=l k=l 

Here Wak{G) takes values Wak{C) = dak^bak^Caki depending on the configura- 
tion C. Because of (2.1), Zn = Zj\f{Xi, . . . , Ajv; z^i, . . . , vn) where A's and v^s are 
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regarded as 'variables'; rj is regarded as a parameter (having the meaning of a 'cou- 
phng constant') and it is often omitted in notations. In QISM the dependence on A's 
and i/'s play an important role (in particular, is invariant under permutations 
within each set of variables). 

2.2. QISM formulation. We now define the main objects of QISM in rela- 
tion to the model. First, let us consider vector space and denote its basis vectors 
as the spin-up and spin-down states 

iT)=(j), u)^(;). 

To each lattice row and column we associate vector space C^. We also use the 
convention that upward and right arrows correspond to the 'spin up' state while 
downward and left arrows correspond to the 'spin down' state. 

Next, let us introduce the quantum L-operator, which can be defined as a 
matrix of the Boltzmann weights. Namely, to each vertex being intersection of 
the a-th vertical line (column) and the fc-th horizontal line (row) we associate the 
operator Lak{^a, v^) which acts in the direct product of two vector spaces C^: in 
the 'horizontal' space Tik = (associated with the fc-th row) and in the 'vertical' 
space Va = (associated with the a-th column). We regard arrow states on the 
top and right edges of the vertex as 'in' indices of the L-operator while those on 
the bottom and left edges as 'out' ones. Explicitly, the L-operator reads 

Lak{K, t'fc) = sin(A„ - ffc + 'nTlal) + sin(2r?)(T-fTj^ -|- t+u^). 

Here r's (cr's) are Pauli matrices of the corresponding vertical (horizontal) vector 
spaces. 

Further, we introduce the monodromy matrix, which is an ordered product 
of L-operators. We define monodromy matrix here as a product of L-operators 
along a column, regarding the corresponding vertical space Va as an 'auxiliary' 
space. The horizontal spaces Tik will be regarded as 'quantum' spaces; the space 
TL = Hk is therefore the total quantum space. In defining the monodromy 

matrix it is convenient to think of L-opcrator as acting u\Va®'H and, moreover, 
writing it as 2-by-2 matrix in Va, with the entries being quantum operators (acting 
in W), 

i«.(A«,.,)=f^i^(^?7/^V^^ sin(2,)a, \ _ 
V sin(2r/)CT+ sm(Ac« -ffe -7?a|)yj^^j 

Here the subscript indicates that this is a matrix in Va and ctj, {I = +,—,z) denote 
quantum operators in Ji acting as Pauli matrices in T-Lk and identically elsewhere. 
The monodromy matrix is defined as 

Ta{}^a) = LaNi^a, Vn) ' ' ' Lq,2(Aq,, V2)La\{^a, ^l) 

_ (A{\a) B{\a)\ 

\C{Xa) D{Xa))^^^^- 

The operators A{X) = A{X; vi,. . . , vn), etc, act in H. They play an important role 
in QISM. 

Operators A{X), B{X), C(A), and L'(A), admit simple graphical interpretation 
as columns of the lattice, with top and bottom arrows fixed. Let us introduce 'all 
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spins down' and 'all spins up' states 

N N 

m = ® iT)fe, i^> = ® (2.4) 

k=l fc=l 

where |t)fc and \i)k are basis vectors of Hk- In the case of domain wall boundary 
conditions each column corresponds to an operator B{Xa) (where a is the number 
of the column) while vectors (2.4) describe states on the right and left boundaries; 
the partition function reads: 

Zn = ■ ■ ■ B{X2)B{X,)\it). 



2.3. Izergin-Korepin formula. In [1] Korepin established recursion rela- 
tions for the partition func:tion, and in [2] Izcrgin showed that these relations are 
satisfied by the following explicit expression (see also [3] for details) 

1 ll<a</3<Ar "■\'^0' "^a) lll<j<k<N "'V^3 ' ^k) 

where M. is N-hj-N matrix with entries 

Mak = V>{^a,yk), ip{\v)=— -, (2.6) 

a(A, v)b{\, u) 

while a{\,v), b{\,v) and c are defined in (2.2), and function (i(A, A'), standing in 
the pre- factor of (2.5), is 

rf(A,A') := sin(A-A'). 

In the next Section we sketch a proof of (2.5), originally given in [12], which uses 

exclusively the Yang-Baxter algebra. The method is useful since it can be gener- 
alized to the case of correlation functions (in contrast to the original approach of 
Korepin and Izergin, [3]). 

In the homogenous limit, i.e., when Aq, = A and Vk = 0, expression (2.5) 
simplifies to 

[sin(A-r?)sin(A + r7)]^' 



where N-hy-N matrix J\f has entries 

sin(2?7) 



detM (2.7) 



J^ak = V(A), ^(A) := v'(A, 0) 



sin(A — rj) sin(A -I- 77) 

Expression (2.7) was given for the first time in [2]; the derivation of (2.7) from (2.5) 
was explained in detail in [3]. 

2.4. Emptiness formation probability. Let us denote by i^^'^^ the prob- 
ability of having all arrows on s first horizontal edges (counted, as usual, from the 
top of the lattice) between r-th and (r-l-l)-th columns, to be all pointing left. Using 
operator formalism we can define this probability as 

p{^,s) ^ Z-^{ii.\B{XN) ■ ■ ■ B{Xr+l) TTl • • • TT, B{Xr) ■ ■ ■ -B(Al)|^). (2.8) 

Here ttj denotes the projector on the spin-down state (which correspondingly fixes 
the arrow to be pointing left). 
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We shall call correlation function (2.8) as emptiness formation probability (EFP), 
adopting the name of similar object from the quantum spin chain context [18]. 

In this respect we comment that the name 'emptiness formation probability' 
for the quantity defined by formula (2.8) has to understood with some care, since 
it is in fact corresponds to a different object. Indeed, if one regards, following the 
common practice in QISM, the vector |-||) as an 'empty' state and B's as 'creation' 
operators over this state, then our definition actually corresponds rather to some 
'fullness formation probability'. Nevertheless, we shall follow the tradition and call 
this quantity EFP. 

One can consider, instead of EFP defined by (2.8), the true emptiness formation 
probability, defining it by replacing all tt's in (2.8) by Wj = 1 — Wj. However, 
contrarily to the standard situation in quantum spin chains in absence of external 
field, such quantity cannot be related to EFP defined by (2.8). This is due to the 
fact that in the domain-wall six- vertex model the local polarization is non- vanishing 
almost everywhere over the lattice, or, in other words, the spin-reversal symmetry 
is broken by the boundary conditions. 

Our clioicx; of correlation function (2.8) is motivated by its further application to 
study limit shapes of the model. Indeed, because of peculiarity of both the domain- 
wall boundary conditions and the six- vertex model rule of conservation of incoming 
and outgoing arrows through each lattice vertex, EFP (2.8) actually measures the 
probability that all vertices in the top-left {N — r) x s sublattice have the same 
configuration of arrows, namely, all arrows point to the left or downwards. Hence 
EFP measures ferroelectric order or 'freezing' of states. The limit shape arises in 
some appropriate scaling limit and corresponds to some curve where EFP jumps 
from one to zero as the size of this {N — r) x s sublattice increases (see [15] for 
further details). 

3. QISM and recurrence relations 

3.1. Yang-BeLxter algebra. One of the most basic relations of QISM is the 
so-called "RLL" relation [18,23,24], which reads 

-Raa'(A, A') [i„fc(A, ly) ® L„.fe(A', l^)] = [Lafe(A', l^) «) La'fc(A, 1^)] Raa'{\ A'). 

Here if!aa/(A, A'), called the _R- matrix, is a matrix acting in the direct product of 
two auxiliary vector spaces, Va ®Va'^ and it can be conveniently represented as a 
4-by-4 matrix (we assume that the first space refers to the 2-by-2 blocks, while the 
second one to the entries in the blocks): 



Raa' (A, A') 



//(A', A) \ 

.9(A',A) 1 

1 5(A',A) 

Here the functions /(A', A) and g{\' , A) are 



/(y,A) = ^^^(.\-/' + ^^) , ,(A',A)= ^^"(2^) 



sin(A-A') ' ' ' sin(A-A') 

This i?-matrix is also sometimes referred to as XXZ chain i?-matrix, due to relation 
of the six-vertex model with Heisenberg XXZ quantum spin chain. It is to be 
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mentioned that here and below we are mainly following notations and conventions 
of book [18]. 

The importance of the RLL relation above is that it implies the following rela- 
tion, which, in turn, can be called RTT relation, 

i?„a'(A, A')[T„(A) ®T„,(A')] = [T„(A') ®T„,(A)]i?„„'(A,A'). (3.1) 

This relation contains all commutation relations between the operators ^(A), B{X), 
C(A), and D{X). The algebra of these operators is called the Yang-Baxter algebra, 
or the quantum algebra of monodromy matrix. Among the sixteen relations con- 
tained in (3.1) the following two of them will be used explicitly below, namely, 

B{X)B{X') = B{X')B{X). (3.2) 

and 

A{X) B{X') = /(A, A') B(A') ^(A) + g{X', A) B{X) A{X'). (3.3) 

3.2. The 'two-site model'. Let us consider the following decomposition of 
the monodromy matrix 

r(A) = T2(A)Ti(A), 

where T'i(A) is defined as a product of the first several L-opcrators while T2(A) 
is the product of the remaining ones. Such a decomposition is sometimes called 
'two-site model' [18]. We shall consider here the case when Ti(A) consists of just 
one L-operator, 

T2(A) =Laiv(A,z/jv)---ia2(A,i/2), Ti(A) =L„i(A,!/i). (3.4) 

Introducing the operators Ai(A), A2{X), etc, as operator- valued entries of the corre- 
sponding monodromy matrices T'i(A), T2(A), respectively, taking into account that 
B{X) = A2{X)Bi{X) + B2(A)£>i(A), and using (2.3), we have 

B{X) = A2(A) c(7f + B2{X) [a(A, !/i)7ri + b{X, vi)Tti] 

^ (h{\ui)B2{X) \ 

V c^2(A) a{\v,)B2{X))^^^^- 

The lower-triangle structure of operator B{X) as a matrix in Wi leads to the 
property that the product of several S's has the form 

B(A.,...B(..) = (!;<;;;;;;:;« ^,(..^^, (3.=) 

The diagonal entries £^ii(Ai, . . . , A„), E22{X\, . . . , A^) are simply proportional to 
B2{Xn) ■ ■ ■ B2{Xi), while the non-diagonal entry, E2i{Xi, . . . , A„), reads 

n ri, f.v — 1 

-B2i(Ai,...,A„) = ^ Yl a(A/3,z^i) • c- 6(A0,j/i) 

a=l 0=a+l 13=1 

X B2{Xn) ■ ■ ■ S2(A„+i)^2(A„)S2(A„_i) • • • B2(Ai). (3.6) 

It is to be stressed that due to (3.2) entries in (3.5) are totally symmetric under 

permutations of A's. While this is completely evident for the diagonal entries, 
such a property is rather non- trivial for expression (3.6), and is a consequence of 
commutation relation (3.3) for operators .62 (A) and ^2 (A). 
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3.3. The key relation. In dealing witli tiie 'two-site model' it is useful to 
consider the corresponding decomposition of the vectors 'all spins up' and 'all spins 
down', e.g., 1^) = |^i)®|^2). To fit (3.4) we set = |T)i and |^2) = <=2lT)fc- 
Obviously, we have 

(^i|S(A„) • ■■BiX.M) = E2i{\i, XnM2). (3.7) 

Taking into account that 

N 

A2(A)|^2) = n«(^''^fc)l^2) 

fc=2 

we can use (3.3) to reduce RHS of (3.7) in terms of -B2's only (applied to vector 
|f|-2)). The result reads 

n n n N 

s2i(Ai,...,A„)iir2) = n n /(^"'^/s) n^^^"''"'^) 

a=l /3=1 f3=l k=2 

X B2(A„)---B2(Aa+l)B2(Aa_l)---B2(Al)|^2). (3.8) 

To get a hint about how this formula can be derived, it is sufficient to look 
at the term a = n in (3.6) which is proportional to Ai(\n)B2{\n-\) ■ ■ • 
This is the only term which contributes to the a = n term of (3.8) (containing 
i3(A„_i) • • •i3(Ai)) after applying commutation relation (3.3); moreover, only the 
first term in RHS of (3.3) contributes to the a = n term in (3.8). The remaining 
terms in (3.8) are just due to the total symmetry with respect to permutations of 
A's. 

Formulae (3.7) and (3.8) express a vector containing B's in terms of vectors 
containing S2's, so they can be seen as a recurrence relation with respect to TV, 
the number of lattice sites. Choosing a specific value of n (the number of B's) and 
specifying also the component of this vector, one can obtain recurrence relations 
for some important 'scalar' quantities, such as partition function and EFP of the 
domain-wall six-vertex model. In turn, these recurrence relations can be explicitly 
solved. 

3.4. Recurrence relation for Zjv. As an illustration of the approach let us 
consider here how the recurrence relation for Z^r can be derived and solved. The 

recurrence relation emerges when n is specified to the value n ~ N, and taking the 
scalar product of (3.7) with the vector (JJ.2|. In this case relation (3.8) gives 

N N N N 

^^ = ^J2]1 H^P'^i) n n a{^a,l^k)ZN-l[Xa;i^l]- (3.9) 

a=l /3=1 /3=1 fe=2 

Here Zjv-i[Aa; i/i] denotes the partition function of the domain- wall six- vertex 
model on (iV — 1) x (A^ — 1) lattice, with the sets of A's and i^'s such that they do 
not contain A^ and vi, namely, they are Ai, . . . , A^-i, Aa+i, . . . , Aat, and ^2, . . . , vn, 
respectively (in other words, the square brackets indicate independence of these 
variables, in comparison with the 'original' sets Ai, . . . , A^r and vi, . . . , vn)- 

Relation (3.9) obviously represents a recurrence relation for the partition func- 
tion with respect to the size of the lattice. The initial condition to the recurrence 
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is Z\ — c. It is to be emphasized that in (3.9) values of A's and z/'s are completely 
arbitrary. 

In order to prove that Izergin-Korepin formula indeed solves (3.9), one can just 
substitute (2.5) into both sides of this relation, and verify whether this is an identity 
or not. After substituting (2.5) into (3.9), and cancelling many factors, one is left 
with 

detA^ = ,nL2^^i^^(_l)-i^(A„)detA^[„;i]. (3.10) 

Here A4[a;i] denotes the {N — 1) x {N — 1) matrix obtained from matrix JVi (2.6) 
by removing a-th row and the first column. The function g{X) is defined by 

^(A) := ^.^'Ifr'l (3-11) 

where 

e(A,A') :=sin(A- A' + 2r/). (3.12) 
It is useful to note that the sum in (3.10) can be written as the determinant 
g{Xi) (/?(Ai,i/2) ... <f{Xi,PN) 



fif(Ajv) ip{XN,1^2) ■■■ ip{XN,l^N) 

The main point is that for each value A = Aa (a = 1, . . . , N) the function g(A) = 
^(A; Ai, . . . , A]v; I'l, ■ ■ ■ , vn), where Ai, . . . , Ajv and v^s are to be regarded as param- 
eters, can be represented as follows 

g{X^) = ^^(A„,.,)i^^^i^kW^. (3.13) 

Noting that the ratio of the two products here does not depend on a = 1, . . . ,N , 
it is easy to see that relation (3.10) is indeed fulfilled. 

The only nontrivial point in this derivation is identity (3.13) which can be 
proven by various methods (e.g., by induction). Let us mention here a hint (even if 
not completely rigourous) on how such kind of relations can be deduced: at values 
A = t'fc + 7? (fc = 1, . . . , A'') function g{X) has simple poles; the sum in RHS of (3.13) 
is nothing but the sum over these poles, in analogy with meromorphic function 
expansion. 



4. Calculation of EFP 
4.1. Recurrence relation. Let us denote 

that is, we will consider temporarily the sole matrix clement of (2.8), or the 'nu- 
merator' of the correlation function. Choosing in (3.7) n = r and taking the scalar 
product of this vector with (JJ-2|52(Ajv) • • •S2(Ar-i-i), one can easily see that (3.8) 
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implies the following recurrence relation 

N r r r N 

fn'^ = n '^(^z'' ^ E n n n ^fc) 

0=r+l c«=l /3=1 /3=1 fe=2 

xF^^ri'^-'^[A„;<^i]. (4.1) 

As above, the 'dependence' on some A„ and Vk enclosed in brackets indicates that 

the sets of A's and z^'s on which the quantity depends, has N — 1 elements each, with 
Xa and missing in the corresponding set (i.e., brackets indicate no dependence 
on these variables). 

In the case s = 1 EFP describes the 'boundary' polarization, which was con- 
sidered in [12]. For all r = 1, . . . , A^, we have 

i^jl'ri'°^[Aa; i^i] = Zj^_,[X^;p,] (4.2) 

and therefore we can just plug Izergin-Korepin formula into RHS of (4.1) in order 
to obtain EFP at s = 1. Indeed, taking into account that 



TT (ll^Ag, /Kg) -p-r 

Zn a{Xa,i^i)b{Xa,i'i) a{Xp,iyi)b{Xp,i'i) ^J-^a{Xa,i^k)b{^a,J^k) 



/3= 



X — 4.3 
detAl ^ ^ 

we arrive at the expression 
p{r,i) _ 1 nf=2 d{iyi,iyk) ex ^ Hpt A/f, , 

aetJVl na=ia(Aa,I/i)na=r+lHAa,I^l) „=! 

(4.4) 

Here the function gr{X) := g(r(A; Ai, . . . , Ajv; J^i, • . • , i^n), is given by 

nl.+i^(Aa,A)nug(Aa,A) 
rifc^i ^(A, '^k) 

where function e(A, A') is defined in (3.12). For r = N function gr{X) is just function 
5(A) defined in (3.11). 

4.2. The cases s = 2 and s = 3. Using (4.1) and (4.4) we can derive EFP in 
the case of s = 2. Indeed, for a = 1, . . . , r, using (4.4), we can write 

AT r N 

detM[„,^ j^J-^ j-ii«(A/3,i^2)^4-^^6(A/3,i/2) 

X j:{-lf^'+^^''"^^i^^,9r{X,) detAl[„,,;i,,j. (4.5) 

Here x(/3, a) = 1 if /3 > a, and x(/3, a) = otherwise. Substituting in (4.5) the 
expression for detA^j^.i] which follows from relation (4.3), and switching to the 
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non-normalized quantity, we get 

;(r-i,i)p _ , _ Zn yT d,{Xf3, Act) TT- d{vi,Vk) 



N r N 

xn'^(^2,^^fc)n^^^ — ^ n — ^ 



e(Aa,A;3) 

Plugging this expression into the recurrence relation, after all cancellations, we 
finally find 

pir,2) ^ n^=2 d{vi,Vk) nr=3 d{v2, Vk) 

" detA< ■ nLia(A„,z^i)a(A„,i/2)nl,.+iKAa,z^i)6(A„,i/2) 
X E ^(-l)"+^+^(^-") '^^^"'"^^^i^^^'"^) g.(A.K(A,)det.Mi.,,^,,,;. (4.6) 

a=l,3=l e(A„,A;3) 

Expression (4.6) can be further used to find EFP for s = 3, by repeating the 
procedure just explained. We quote here only the result. 



A< fA nUl i',) nl^+l KAa, J^,) 



■'V A/< 11 T-rr ,M-rJV 

i 

r r 

X 



E E E (-l)"+''+"+'+^''''"^+^'"'''^+^^^'"^5.(A«)ff.(A0)<?.(A,) 

- 7=1 

a(Aa, z^2)a(Aa, j^3)a(A/3, j^3)&(A;3, :/i)&(A^, i/i)6(A^, 1^2) 



a=l /3=1 7=1 



e(Aa, A/3)e(Aa, AT,)e(A/3, A^,) 

X detA1[„,/3,^.i,2,3] . (4.7) 

4.3. Result for generic ,s. Inspecting formulae (4.4), (4.6) and (4.7) it is 
rather straightforward to guess the result for generic values of s. The following 
expression for EFP is valid: 



1 -Q Y{k=j+ld{^j^^k) 



'n ^ det li n^^i a(A^, u,) HJ^.+i 6(A^, v,) 

r r r s 

X E E ••• E (-i)^+^^=i"'=+2:i<^<'=<='^("'='"^) [|5^(A„j 

ai — la2— 1 cts — 1 i— 1 

^ TT «(Aa,,t/fc)6(Aa^,t/j) 

X 11 TT detA1[„j,...,„^.i,...,s] . (4.8) 

It can be shown directly that formula (4.8) solves recurrence relation (4.1). The 
calculation in fact repeats the one above for the case of s = 2. 
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Let US focus on writing EFP for the case of (A'' — 1) x {N — 1) lattice, namely, 
we are interested in -F^ri '''~^^[A„^, i^i]. To be more precise, in adapting expression 
(4.8) to this case, apart from using the fact that the sum is now (s — l)-fold, it is 
important also to take into account that the function gr{^) implicitly depends on 
N, so that switching from N to {N — 1) implies some extra factor. In all, after all 
these preparations, which arc mostly devoted to fit the notations, the expression 
which has to be substituted into RHS of (4.1) reads: 



^ir-l,s-l) 
JV-1 



[Aai , Z^l] = 



n 



det A^[„,;i] 11 n'^^i «(A/3, Vi) YV;=r+i HA/3, Vi) 



^ ... ^ (_1)«-1+Efc=2«fc+Ei<3<ji<, X(afe,aj) 



^2 — 1 CVs — 1 

a2 7^0ii Q,y7^^Qi,...,a 



\ ai, ajj 2<j<k<s 



a(Aa,,i/fc)b(AQ^,t/j) 
e(Aa,, AaJ 

X det7V4[„^,...,„^.i,, 



(4.9) 



Again, as in the case of s = 2 considered above, one can use relation (4.3) to elim- 
inate det in favour of det in (4.9). With this point taken into account, 
cancellation of various factors shows that recurrence relation (4.1) is indeed fulfilled. 

Apart from showing that (4.8) is a solution to (4.1), it is also useful to mention 
that the corresponding 'initial conditions' are also satisfied. Indeed, for s <r, such 
a condition is just the generalization of condition (4.2), namely 

F^-s ' ''[Aai , • • • , Aas ',l^l,---,l^s\ = ZN-s[Xai , • • • , Aqis ', 1^1, ■ ■ ■ , l^s]- 



For s > r, the corresponding condition is 



F^-r ^ [-^an • • • ! ^ctr ; , • • • 5 i^r] — 0. 

Note, that this relation implies the relation F^'*-* = which must hold whenever 
s > r, and which follows directly from the definition of EFP (2.8). 



5. EFP in the homogeneous limit 

5.1. The procedure. The homogeneous limit can be performed along the 
lines of papers [3, 14]. We start with writing 

Aa = A + , 

where ^'s will be set equal to zero in the limit (as well as i^'s). Keeping ^'s nonzero 
(and different from each other), and using the fact that for a function f{x), regular 
near x = X, one has exp(^9£)/(A-|-£)|e=o = /(A + ^), we can bring (4.8) into a form 
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which involves some determinant: 



[7'(''.«) 
■ JV 



1 



n 



exp(^i(9ej 



xl[gr{X + ej) n 

j=l l<j<fe<s 



exp{S,Nde,) tpiXN, Ks+i) ■■■ v{Xn,'i^n) 
a{X + ej,Vk)h{X + e^, i/j) 



(5.1) 



e(A + ej,A + efe) 

£i=...=e,,=0 

It is to be emphasised that this expression is vahd for the inhomogeneous model 
(no homogeneous limit yet); it represents an equivalent way of writing the multiple 
sum in (4.8). 



Let us now perform the limit v\ ^0, V2 



0, 



0, in this order, at 



each stage keeping the contribution of leading order in the corresponding variable; 

next wc do the same with ^'s, in the order 0, ^2 ^ 0, . . ., 0, again 

keeping only the contributions of leading order. The first line of (5.1) gives 



(-1) ' Uj^A^-jV- 

a"6(^-'-)MetA/' 



1 



1 



1! 



2! 



(Af-s-1)! 1! 2! 

where a := a(A,0), b := 6(A,0). The second line of (5.1) gives 



£2. . ia. . . . ''w 



(JV-l)! 



(-1) 



(N- 



-l)(N-s) 



1 



1 



(^(A) 
dx^iX) 



d 



^(A) 
9--V(A) 



JV-l aN-l 



ar-v(A) 



,N-s-l 



\ 1! 2! (iV-s 
The terms in the third line of (5.1) give 



1)! 



ar-^-v(A) 



6 ei 

1! 2! (A^-l)! 



(5.2) 



-Q [_- sin(£j- - 2ri)Y{- sinej)^" 



[sin(£j + A - 77)] 



N 



n 

l<j<fc<s 



sin(£j + A + T/) sin(e/j; + A - 77) 
sin(ej - efe + 2?7) 



In all, after all cancellations, EFP in the homogeneous model reads 



■ JV 



(-irn;=i(^-j)! 

ars5(JV-r)sdet7V 



^(A) 
dxifiX) 



d 



N-s-l 



^(A) 



1 



2N-S-2 



5r-V(A) ... a; 

(sin£,)^-'-[sin(e,-2r7)]^ 



JV-l 



1 



[sin{ej + A - 77)] 



N 



X n 

l<J<fe<S 



sin(£j + A + 77) sin(efc + A — 77) 



(5.3) 



sin(e,- - Efe + 277) 

£l—...—S,^—0 

Note, that in the determinant here we have changed the order of columns, in com- 
parison with formulae (5.1) or (5.2). 



EMPTINESS FORMATION PROBABILITY 



13 



5.2. Orthogonal polynomials representation. Formula (5.3) for EFP of 
the homogenous model, involving N x N determinant, can be also represented in 
terms of some s x s determinant. Such an equivalent representation can be further 
used to obtain a multiple integral representation (having s integrations), similar to 
those for correlation functions of quantum spin chains [19-22]. 

The derivation of the s x s determinant representation from (5.3) is based on 
the following facts. Let {Pn{x)}^^o be a set of orthogonal polynomials. 



/ 



Pn{x)Pm{x)lJ.{x)dx = hn6„ 



(5.4) 



where the integration domain is assumed over the real axis, the weight iJ,{x) is 
real nonnegative, and we choose ft.„'s such that Pn(a:)'s are monic (i.e., the leading 
coefRcient is equal to one, Pn{x) = + . . . ). Let c„ denote n-th moment of the 
weight n{x), 

x"/i(.T)d,T (n = 0, 1,...). 

The orthogonality condition (5.4) and standard properties of determinants allow 
one to prove that 



Co Ci 
Cl C2 



Cn-1 



hghi • ■ ■ hn-1- 



Cn-1 C„ ... C2ra-2 

More generally, the following formula is valid 



Co 


Cl . 


Cn—k—1 


1 


1 


1 


Cl 


C2 • 


Cn—k 


Xl 


X2 


Xk 


Cn-1 


Cn 


■ C2n-k-l 


xr' 


X2 





ho - ■ ■ hn-k-l 



P„-k{xi) ■ 


■ Pn-k{xk) 


Pn-lixi) . 


■ Pn-l{Xk) 



(5.5) 



The proof of this formula can be given by the methods explained in [25]. We shall 
use (5.5) with n = N and k = s. 

To apply (5.5) to expression (5.3) we need first to identify the orthogonality 

weight. Since in our case, by definition, c„ = 9"(p(A), the weight can be obtained 
by representing the function (p(X) via Laplace transform; 



V(A) 



e^^<&(a;) da;. 



Here function $(a;) = ^{x;t}) is independent of A, but depends on 77 as a param- 
eter (indeed, in our case /i(.T) = fi{x;\rj) where A and 77 are to be considered as 
parameters, and ijl{x) = e'^^<^{x)). 

The explicit form of function $(x) depends on the particular choice of the 
physical regime in the six-vertex model. For example, let us consider the so-called 
disordered regime (|A| < 1). In this case A and rj in our standard parametrization 
of weights a = sin(A + 77), 6 = sin(A — 77), and c = sin(277) are both real and satisfy 
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< 77 < 7r/2 and 77 < A < tt — 77. Function ^(x) can be found from the formula: 

da; = . r. (5.6) 



r 



.x(A-7r/2) sinh(77 a:) _ sin(277) 



sinh(7rx/2) 



sin(A — 77) sin(A + 77) 



The other two regimes, the ferroelectric (A > 1) and anti-ferroelectric (A < — 1) 
can be approached when A and r] are complex, satisfying Re A = Re 77 = or 
ReA = RcT/ = 7r/2, respectively. In fact, for these regimes the weight ^{x) can be 
found from (5.6) using the proper analytical continuation in the parameters A and 
77. The integral in RHS of (5.6) in these cases is replaced by a sum coming from 
the simple poles of the integrand, or, in other words, the corresponding measure 
^{x)(ix turns out to be discrete (see [5] for explicit formulae). For our construction 
below the actual choice of the regime, and hence the explicit expression for the 
orthogonality weight /i(x), is irrelevant since only the corresponding orthogonal 
polynomials will enter the formulae explicitly (for instance, one can think of the 
case of the disordered regime; the resulting expressions in terms of P„(a;)'s for other 
regimes turn out to be essentially the same). 

To make the use of (5.5) convenient in the framework of representation (5.3), 
we define, following [14], the functions 



He) 



sin(A + 77) sin e 



sin(A — 77) sin(£ — 277) ' 



a;(e) 



sin(A — 77) sine 



sin(A + 77) sin(e + 277) " 



(5.7) 



They are functions of e, with A and 77 regarded as parameters, as indicated in the 
notations (below we omit the dependence on A and 77 whenever possible). We also 
define 



sin(A — 77) sin(e — 2r;) 



sin(277) sin(e + A — 77) ' 
There are useful relations 

1 



u;{e) 



sin(A + 77) sin(e + 2r;) 



sin(277) sin(£ + A + 77) " 



l-u;{e) 



(5.5 



Noting that 



siu(ji + A + ;/) siu(j2 + A - 77) 



and defining 



sin(£i - £2 + 277) 

Kn{x) = 



ipQ{£l)Q{e2) Cj{£i)ljj{e2) - 1 
Pn{x) 



(5.9) 



we have, in virtue (5.5), the following orthogonal polynomials representation: 



■ N 



X n 



l<j<k<s 



£i=...=es=0 



(5.10) 
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5.3. Multiple integral representation. Let i?]^ denotes the probability 
of having the c-weight vertex of the first row (from the top of the lattiee) at r-th 
position (from the right). This quantity can be related to EFP at .s = 1. wliicli 
describes in this case essentially the one-point correlation function (polarization) 
at the boundary, via H^^ = F)^'^^ — F^~^'^^ (we refer for details to [12] where 

Fj^'^^ was denoted as G^^^). Using (5.10) and taking into account the first relation 
in (5.8), we have (see also [14]) 



(5.11) 

e=0 



Let us define the generating function: 



N 



r=l 

The key identity, for the derivation of a multiple integral representation for EFP, 
is 

KN-x{d,) f{u^{e)) ^-^^^^^hN{z)f{z)<iz. (5.13) 



e=o 27ri 



Co 



Here f{z) is some function, regular near z = Q, and Co is a simple closed counter- 
clockwise countour around the origin. 

To prove this identity, we note first that since J{z) is assumed regular at z = Q, 
one can further think of it as a polynomial of degree {N — 1), since higher powers 
in z do not contribute to either side of (5.13) (note that w{e) — > as £ — > 0). It 
therefore suffices to prove (5.13) for f{z) being a monomial. Next, let us define 
auxiliary quantities 

y^^) = i^jv-i(9e) Her-\^^, (p = 1, . . . , TV). (5.14) 

Evaluating the integral in RHS of (5.13), one finds that (5.13) is nothing but the 
relation 

^={-l)''-^Ah, (5.15) 

where the vectors v and h have components Vp = Vjif' and hr = H)^' , and matrix 
A has entries Apr = {^^Y~^{plr)- Obviously, matrix A is lower-triangular, and 
moreover it can be represented as 

A={I-£)^-\ (5.16) 

where £ denotes the lower triangular matrix, with entries standing under the main 
diagonal equal to one and all other entries being zeroes, i.e., £pr — Sp_r+i- Inverting 
matrix A in (5.15) with the help of (5.16), one arrives immediately to the expression 
for H^^ in terms of vjf^^s, which follows from formula (5.11) and definition (5.14). 
This proves identity (5.13). 

As a result, applying identity (5.13) to (5.10) and taking into account relations 
(5.8) and definitions (5.7), we obtain the following multiple integral representation 
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for EFP: 



1 



Co 



Co 



z{{zi - 1) 

/lAf-l(zi) 

zr\z, - If 

hN-s+l{zi) 



ZliZs - 1) 
flN-liZs) 

1)2 



r— 1 / 
Zs {Z. 



r — s+1 



(zi - 1)^ 



r— s+1 



X n 

l<j<k<s 



{Zj - l){Zk - 1) 



ZjZk 



Here Zj's are functions of the corresponding zj^s, 

Pz, 



(a^ + 62 - c'^)zj 



dzi---dzs. (5.17) 



(5.18) 



where a, b, and c are the homogeneous six-vertex model weights (formula (5.18) 
follows from (5.7) by expressing function ui{e) in terms of w(£) and using a = 
sin(A + 7]), b = sin(A — jj), and c = sin(2ry)). 

6. More on the multiple integral representation 

6.1. Preliminaries. Formula (5.17) gives EFP in terms of a multiple integral. 

This integral representation can be transformed into some other forms which can 
be useful for further study of EFP (e.g., for finding the limit shape). Our aim in 
this section is to derive some of these representations. 

The first point to be mentioned in this respect is that the integrand of (5.17) 
involves some determinant, which is an antisymmetric function with respect to 
permutations of the integration variables Zi, . . . ,Zs- Because of this antisymmetry, 
only the antisymmetric part of the double product, with respect to permutations 
of these variables, actually contributes to the multiple integral in (5.17). 

Let us introduce the parametrization 



Zj = ^i^j), Zj = Qi^j) 



(6.1) 



where the functions w(^) = uj{^;X,ri) and a)(^) = a;(^;A,7?) are as in (5.7) (note 
that upon this parametrization relation (5.18) is satisfied automatically). With 
this parametrization the double product in (5.17) coincides with the expression 
whose antisymmetric part has been found in paper [20] (see appendix C of that 
paper). Namely, using a mixed set of notations, the antisymmetric part of the 
double product, with respect to permutations of the variables zi,...,Zs, can be 
written as 



Asym JJ 



" l<j<k<s 



{Zj - l){Zk - 1) 
ZjZk - 1 



i<s</'' " .4 ^^-^-^ - - + 



EMPTINESS FORMATION PROBABILITY 



17 



Here Zs{Xi, . . . , A^) denotes the partition function of the 'partially' inhomogeneous 
six-vertex model with domain wall boundary conditions on the s x s lattice; by 
'partially' inhomogeneous we mean that the weights of this model are given by 
(2.1) and (2.2), but with all f's equal zero. 

The expression in the last line of (6.2) is some symmetric polynomial in variables 
zi,. . . ,Zs. Our aim below will be to show that this expression admits a representa- 
tion in terms of some determinant. Furthermore, due to this representation it turns 
out possible to clarify the meaning of the determinant in (5.17) which appears to 
be related to the partition function of the partially inhomogeneous model on the 
N X N lattice with s inhomogeneities (out of N possible). 

6.2. Back to the partition function. Let us consider the 'partially' inho- 
mogeneous six-vertex model on N x N lattice, whose partition function will be 
denoted as Zn{Xi, . . . , Xn)- The weights of the model are given by (2.1) and (2.2), 
with fi = 1/2 = ■ ■ ■ = vn — 0, while Ai, A2, . . . , Ajv are kept in general to be differ- 
ent. Choosing some particular A as a value for Ai, A2, • . • , Xn in the homogenous 
limit, one can introduce variables ^1, ^2, • • • , ^jv by 

A« = A + ^„, (6.3) 

i.e., such that the homogenous limit will correspond to making all ^'s vanish. To 

simplify notations wc shall write a{Xa) := a(Aci,0) and h{Xa) '■= 6(Aq,0). Since 
both the homogeneous and inhomogeneous models will be considered here simul- 
taneously, we shall write a(A) and 6(A) for the weights of the homogenous model 
(where A is the same as in (6.3)), rather than simply a and b. We come back 
however to these simplified notations in the last subsection. 

An important role below will be played by the variables ui, . . . ,un, defined by 

For later use wc mention here that function 7(^) is related to function u){^) (5.7) 

via7(0 =w(-A-^^?-0• 
Our aim here will be to study the quantity Z^i^i, ■ ■ ■ , ^n), which depends also 
on T] and A as parameters, and is defined by 

Zn{Xi, . . . , Ajv) = Zn{X, . . . , A)Zjv(Ci, . • . ,^n)- 

As we shall see, there is a nice representation for this "bare" partition function in 
terms of variables (6.4). 

To illustrate the idea, it is useful to consider first the case when only one 
inhomogeneity is present. Without lack of generality we can assume that this is the 
inhomogeneity of the first row (i.e., ^2 = • • • = ^at = 0). The domain wall boundary 
conditions admit only one vertex of weight c in the first row; if this vertex is at 
r-th position (counted, as usual, from the right) then the first (r — 1) vertices are 
of weight b while the remaining {N — r) vertices are of weight a. Taking this into 
account one has therefore the expression 



^jv(Ai,A,...,A) = ^ 

r=l 



N r 



r«(Ai)i 


N-r 


p(Ai)l 


[ a(A) J 




[ 6(A) J 



H)^>{X,...,X)Zn{X,...,X). 
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Note that here hI^\x, . . . , A) denotes the correlation function of the homogeneous 
model (recall that this function describes the probability of having the c- weight ver- 
tex at r-th position on the first row). Recalling definition (5.12), this last expression 
implies that 



zw(a,o,...,o) 



a(Ai) 



a{X) 



N-l 



/lAr(ui). 



(6.5) 



It is to be emphasized that the variables ^i, Ai, and ui are related to each other 
by (6.3) and (6.4); A is to be regarded as a parameter, entering also function 7(^) 
in (6.4). 



To consider the general case, let us introduce functions liN.siui 



where 



the second subscript in the notation refers to the number of arguments (s 
1,...,N) 



hN,s{ui,...,Us) = {Uk-Uj) ^ 

l<j<k<s 

ul~'^hN-s+l{Ul) 
Wl~^(wi - l)/ljV-s+2(ui) 

(«i - iy-^hN{ui) 



ul~'^{Us - l)hN-s+2{Us) 



(6.6) 



These functions are symmetric polynomials of degree (TV — 1) in each of their vari- 
ables. They satisfy the relation (recall that h^i^) = 1) 



hN,s+i{ui, • • • , Ws, 1) = hN,s{ui, ■■ ■,Us). 



(6.7) 



This relation says that these functions can be in fact constructed iteratively, starting 
with /ijv,jv(wi, . . . , ujv); one has then, in particular, hN,i{ui) — /lAr(ui). The "bare" 
partition function, in the most general case of N inhomogeneities, reads 



Zn{^i, ■ ■ • ,Civ) 



N 

n 



a{Xj) 



-I N-l 



a(A) 



hN,N{ui, . . . ,Un)- 



(6.8) 



Using relation (6.7) one can easily obtain the corresponding formula in the case of 
s (out of A'') inhomogeneities. 

It is worth to mention that /ijv,s(mi, • ■ • , w.s) introduced in (6.6), besides describ- 
ing the partition function with s inhomogeneities, turns out to coincide, modulo a 
trivial factor and the substitution Uj i— > Zj, with the determinant in (5.17). We will 
exploit this observation below in discussing equivalent multiple integral represen- 
tations for EFP. 

6.3. Proof of the formula for Zn{^-i, . . . Representation (6.8) for 

Zn{^i, ■ ■ ■ tS,n) (for the special case of A = 7r/2) appeared for the first time in paper 
[26]; since the proof in that paper was lacking, we sketch it here for completeness. 
Let us begin with Izergin-Korepin formula, specialized to the case oivi = ■ ■ ■ = 
=0, which reads 



nf=i[«(A,)&(A,)] 



Zjv(Ai, . . . , Ajv) = 



N 



Hi 



<j<k<N 



N-l . 
n=0 



'/'(Ai) 



(^(Ajv) 



C"V(Aiv) 
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Taking into account that ip{X) = c/a{X)b{X), and passing to the orthogonal poly- 
nomial representation (see formulae (5.5) and (5.9)), we get 



Zn{^i, • • ■,S,n) 



b(A)] 



Af(JV-l)/2 



nf=ik(A,-)]^ ni<i<.<iv diXk, xj) 



Now, taking into account that 

Uk - Uj 

and 



[a{X)f 



Kn -lidxf,) 



N 



a{Xj)a{Xk) 

ip{X)a{X) o,{Xj 



ip{X)d{Xk,Xj) 



(p{Xj)a(Xj) a{X) 



J, J, 



(6.9) 



(6.10) 



we arrive at 



N 



^jv(a,---,^iv)=n 



a{Xj) 



a{X) 
1 



<j<k<N 



{uk - Uj) 



nf=i^(A,) 



Ko{dx,) 
Kiidx,) 



Kiidx.) 



KN-l{^X^) 



KN-i{dx^) 



N 



Yl^iX,). (6.11) 

j=i 



Comparing with (6.8), we see that to complete the proof we only need to express 
the second line in terms of w's. 

For this purpose let us again turn to the special case when only one inhomo- 
geneity is present. In this case Izergin-Korepin formula boils down to 



Zj\f{Xi,X, . . . , A) 



[a(Ai)&(Ai)]^[a(A)5(A)]^(^-i) 



(A^-l)!nt"o'(«!)^M(Ai,A)] 



N-l 



^(A) ... a;-V(A) ^(Ai) 
dxviX) ... 9r-V(A) dxMXi) 

5f-v(A) ■ar-v(A)''<-V(Ai) 



We have 



Zw(ei,0,...,0) = 



^ KN-i{dxi)ip{Xi). 



[(^(Ai)]^KAi,A)]^- 

Taking into account that (6.9) implies that d(Ai,A) = {ui — l)a{Xi)[a{X)(p{X)]~^ 
and also using (6.10), we have 



Ziv(Ci,0,...,0) 



ajXi) 
a{X) 



N-l 



Ml 



Ui — 1 



N-l 



<p(Ai) 



KN-i{dx,MXi). 
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Now, comparing this expression with (6.5), wc find that 

^i^„(a,,)^(A,)=(^) K^,{uj). 

Finally, plugging this formula into (6.11), one arrives at (6.8), which is thus proven. 

6.4. Equivalent representations. Coming back to the expression in the last 
line of formula (6.2), one immediately finds, using (6.8) (and also taking into account 
(6.1), (6.4), and (5.8)) that 



ITn(^^j Ksiui,..., Us). (6.12) 



n;=i[sin(e,-2r/)]-i c^(« 

Here Zg stands for Zs{X,..., A) and below we also come back to our simplified no- 
tations (a = a(A, 0), b = b{X, 0)). Furthermore, we shall use the following notations 

+ 6^ - ,„ , b sinfA - n) 

A:= — -— = cos 2n), t := - = ^4-. r- 

2ab ^ a sin(A + r?) 

Exploiting formulae (5.7) and (6.4) it is straightforward to find that 

Z," — 1 

II - — ^ 

' (^2 - 2tA)zj + 1 ■ 

This formula, by the way, shows that (6.12) is indeed a polynomial in 2:1,..., 2;^ 
(recall that in general the functions h^^siui, ■ ■ ■ ,Us) are symmetric polynomials of 
degree (iV — 1) in each variable). In all, for the antisymmetrical part of the double 
product in (5.17), we obtain 

Asym n ^^^-^^^^ = T#T)T H 

riLiKi' - 2tA)z, + l]^-i hsAui, ...,Us) 

X ^ 

Taking account that the determinant in (5.17) can be represented in terms of the 
function hN,s{zi, ■ ■ ■ , Zs), we arrive therefore at the following equivalent multiple 
integral representation for EFP 



(-i)-(^+i)/2z, / / TT / ^2 

J Co J Co i<j^k<s 

1 A [{t^ - 2tA)zj + l]^-i 



X 



t'^j^k - 2tAZj + 1 n ,r^,. _ 1). 

X hN^s{zi, . . . , Zs)hs^siui, . . . , Ws) dzi • • • dzg. (6.13) 

This formula is the main result of the present paper. 

As a comment to this result, let us also mention that due to (6.8) one can 
obtain from (6.13) a multiple integral representation, in which the integrand is 
expressed in terms of the partition functions. Indeed, since z/s and u/s are given 
by (5.7) and (6.4) where functions u){^) and 7(^) appear to be related to each other 
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as 7(^) = ui{—X + 7] — £,), or ui{£^) = 7(— A + — one may prefer to use ^i, . . . 
as the integration variables, and furthermore (6.8) imphes that 

r f\\ iN-1 

. / /e ^ re \\ ^N\.il-Zi,---,il-';s,^,---,^) 

h^A^{^i),---Mis)) = z^{x,...,x) n 



a(A) 



As a result, we obtain the following alternative multiple integral representation 



X 



S s 



fc=i - + 1=1 - 2r?)]^-1sin(^,- + A - 7?)]nsin^,)'- 

x^jv(??-6,---,^-^.,A,...,A)Z,(A + a,---,A + 6)d6---d6- (6-14) 

This formula evidently recalls multiple integral representations which appear for 

correlation functions in quantum spin chains (sec, e.g., [19-22]). 

The important and essentially nontrivial object in the integrand of (6.14) is the 
partition function Zn, with s shifted arguments. While the presence of the second 
one, Zs, is a standard consequence of the Yang-Baxter algebra, the first partition 
function, Z^, together with other factors in the integrand, is due to the specificities 
of both the boundary conditions (domain-wall ones) and the particular correlation 
function under consideration (EFP). 

In conclusion, let us briefly discuss formula (6.13). This representation may 
seem not completely explicit since it is heavily based on the use of the generating 
function hi^{z). We expect however representation (6.13) to appear most useful 
when addressing the problem of limit shapes in the model. In the case of vertex 
weights obeying the free-fermion condition {q = it/ 4:), function hN{z) is known ex- 
plicitly and particularly simple. This case has been considered in [15]. Fortunately, 
function h^iz) is also known at = 7r/6 and 77 = 7r/3, for A = 7r/2. The first 
case corresponds to enumeration of alternating sign matrices and the second one to 
a particular example of their weighted enumeration (the so-called 3-enumeration). 
For this reason we expect that formula (6.13) may appear also very useful out of 
the free-fermion case, and, at least, allows one to solve the long-standing problem 
of the limit shape of large alternating sign matrices. 
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